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Abstract 

We consider solitonic solutions of coupled scalar systems, whose 
Lagrangian has a potential term (quasi-supersymmetric potential) 
consisting of the square of derivative of a superpotential. The most 
important feature of such a theory is that among soliton masses there 
holds a Ritz-like combination rule (e.g. M12 + M23 = Mis), instead 
of the inequality (M12 + M23 < M13) which is a stability relation 
generally seen in N=2 supersymmetric theory. The promotion from 
N=l to N=2 theory is considered. 
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1 Introduction 



The (l+l)-dimensional solitons are worthwhile problems because the dynam- 
ics of (3+1) -dimensional domain walls are directly connected to them. In 
Ref. [1], attractive intersoliton force is concluded to exist, since there is found 
a mass inequality (M i2 + M 23 < Mi 3 ), where mass refers to a soliton 
interpolating ith and jth potential minimum. It was the N=2 supersymmet- 
ric model which came from D=4, N=l supersymmetric Wess-Zumino model 
through dimensional reduction. In the same model, Ref. [2] and [3] pointed 
out that there is a curve of marginal stability in parameter space where mass 
equality (M 12 + M 2 3 = M 13 ) holds, and absence of the intersoliton force of 
leading order is checked right on the curve. 

Therefore it will be interesting if we have, not as exceptions but as a 
rule, marginally stable bound states of solitons so that the Ritz-like mass 
combination rule [Mij = \Hi — Hj\, where {Hi} is some discrete series) is 
satisfied. Indeed in the cases which we will investigate in this paper, it will 
be a common feature that continuously degenerate class of bound solitons 
comes out. 

As for exact soliton solutions of relativistic nonlinear field equations, ef- 
forts have been continuing since the seventies. Very sophisticated calculations 
are invented to study their properties [4] [5] [6] . Nevertheless, even in two di- 
mensions constructing localized solutions of two or more coupled fields are 
not easy because general systematic method is unknown. The invention of 
supersymmetric theory brought a breakthrough for this purpose [7] . 

We will concentrate on those Lagrangian field theory whose potential 
term is a square of derivative of a certain function, the latter function being 
called superpotential. In the following, we call them quasi-supersymmetric 
as this kind of systems can be made supersymmetric by introducing the same 
number of Dirac fields according to popular prescription. 

In Ref.[l], it was shown that for the D=2, N=2 supersymmetric theory of 
only one superfield, exact soliton solutions are obtainable through Hamilton- 
Jacobi method, and that their masses satisfy a triangular mass inequality. 
This model, actually, is known to be a perturbed conformal field theory. By 
virtue of the solvability, the elastic S-matrix elements were also calculated 
[8]. In contrast to the N=2 theory which is so strongly constrained and has 
fairly known properties, the dynamics of N=l solitons clearly needs different 
approach. In this paper, we investigate such elaborate cases of D=2, N=l su- 
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persymmetry in which two scalar fields couple through quasi-supersymmetric 
potential 



v(4>) 



1 2 



d 

d<pi 



In Sec. 2, starting from (l + l)-dimensional Lagrangian with N=l super- 
symmetry we recapitulate soliton solution. In Sec. 3, we investigate similar 
system with N=2 supersymmetry, characterized by a holomorphic superpo- 
tential. In Sec. 4, we turn to coupled theory of two scalar fields with N=l 
supersymmetry, and derive the mass combination rule. In Sec. 5, the relation 
of N=l and N=2 supersymmetry in scalar soliton theory is considered. 



2 Solitons of Single Scalar Field with N = 1 
Supersymmetry 

We consider the two-dimensional N=l supersymmetric case [9] in which only 
one type of supersymmetric multiplet takes place. The action is 



S = Jd 2 x [-i(^) 2 - \(H'f + jijrfdrf ~ \H"W 



(1) 



where is a real scalar field and if) is a two component Majorana field. H((f>) 
is a real superpotential, H' and H" are its derivative with respect to <fi. 
Equations of motion for <fi and -0 are respectively 

d^<t> - H' ((f)) H" ((p) - \H"\^ = (2) 
i-fd^ - H"((f>)ij = 0. (3) 

Minimum points of the scalar potential V = \(H') 2 coincide with the sta- 
tionary points of H which satisfy equation 

H'(4>) = 0. (4) 

If there are m such points 0(i) > (f)(2) > • ■ • > (f)( m ), they are each candidates 
for asymptotic values 0(±oo). 
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Whenever m > 2, this model has soliton solutions. A soliton is the 
solution of equation of motion which take different values in x — > — oo and 
x — > +00. If these values are 



for 
for 



x 
x 



—00 

+00, 



(5) 



let us call this (i, j)-soliton where i — j — ±1. 

We want to derive the Bogomol'nyi bound of the Hamiltonian. For the 
static solutions, the Hamiltonian is 



n 



dx 



-(H'f 



We can deform this Hamiltonian to 



Ti 



J — i 



dx 



-{d x cj>TH'Y±d x H 



(6) 



(7) 



where we can use either upper sign or lower sign. Since the first term is 
always positive, the lower bound of the Hamiltonian is given by the second 
term. After integrating the second term over x, we obtain the Bogomol'nyi 
bound 



n > ±T (iJ) (8) 

where is the difference of H at (f>^ and (f>^ 

— H{4>{j)) — H((j)(i)). (9) 

|T(jj)| may be called the classical soliton mass M^. On this bound scalar 
field (f) is saturated by solutions of equation 

d x (j>^H' = Q. (10) 

We can confirm that solutions of this equation satisfy equation of motion (2) 
with ip = by differentiating it by x and using it twice. 
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3 Soliton in N = 2 Supersymmetry 

Now we consider the two-dimensional N=2 supersymmetric case [1][2]. The 
action is derived from the four-dimensional Wess-Zumino model via dimen- 
sional reduction. 



s = J d 2 x\-{d^){d^) - (wy(W) 



+i$>fd^ - \w"^ - l -(w"y^% (ii) 

where p is a complex scalar field and ip is a Dirac field. W(<p) is a complex, 
holomorphic superpotential. A soliton is the solution of the equation of 
motion which takes different values for x — > — oo and for x — > +00. This is 
similar to the above section, except now p> is a complex scalar field. Therefore 
minimum points of the scalar potential V = (W')*(W) are distributed on 
the complex plane of p. Solutions of V = coincide with those of 

W'(p) = 0. (12) 

Suppose that this equation has m solutions p^ , p{2) , • ■ ■ , P( m )- When 
i 7^ j boundary conditions of (i, j)-soliton are 

p = </?(j) for x — > —00 

p = if (j) for x — > +00. (13) 

We will derive the Bogomol'nyi mass bound in this model. Equations of 
motion for p> and tp are respectively 

d^p - W'(W"Y - i(VT') W = (14) 
i^d^ip - (W")*^ c = 0. (15) 
When we consider the static solution, the Hamiltonian becomes 

/oo 
dx[(d x p*)(d x p) + (w'y(w% (i6) 
-00 

With arbitrary phase a, we can express this as 



/oo r , 

dx \\d x p - e ia (W')*\ 2 + d x {e~ ia W + e ia W*}] . (17) 
-oo L - 1 



-00 
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Since the first term is always positive, a lower bound of the Hamiltonian is 
given by the second term. For (i, j)-soliton we have 

H>2M[e- ia T {i>j) ], (18) 

where is a topological charge 

— W(if(j)) — W((p({)). (19) 

For this lower bound, scalar field (p is saturated by solutions of equation 

d x (p - e ia (W')* = 0. (20) 

We can confirm solutions of this equation satisfy equation of motion (14) 
with xjj = by differentiating it by x and using its complex conjugate. 
Multiplying this equation by W and integrating over x, we have 

/oo 
dx\Wf. (21) 
-oo 

Thus e ia equal the phase of T^jy If we take 

T(iJ) = W(<pw) - W( m ) = e id \W( m ) - W( m )\, (22) 
the only condition in which solutions of equation (20) become finite is 

e ia = e ie . (23) 

Then mass bound (18) is 

H>2\W( m )-W(<p (t) )\. (24) 

Thus the mass bound is twice the distance of two points on complex W plane. 
Now let us replace 

e - w W = W (25) 

and decompose <p = -^m{4>R + ^4>i) an d W — \{Wr + iWi), where 4>r and Wr 
are real part of f and W, and 0/ and Wj are imaginary parts respectively. 
Further using Cauchy-Riemann equations 

dW R _ dWj_ 
d(p R d$i 
dW R dWj 

dh d<p R [ ' 
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for (20), we have equations 



Ox<PR = -TT— 
(90 



R 

P> A 9Wr (o^ 

o<Pi 

They mean that the integral curve of Eq.(20) coincides with a gradient curve 
of Wn(<pit, (pi). Thus to derive the finite solution, point ip^ must be con- 
nected with <p(j) by a gradient curve of Wr. 

On the other hand, since Wr is a harmonic function, all stationary points 
are Saddle points. We illustrate below that a stationary point is connected 
with other point by a gradient curve of Wr, but it is not connected in general 
by any gradient curves of Wr. 

Let us take the cubic superpotential whose general form is 

W(<p) = cap 3 + b<p 2 + op, (28) 

where a, b, c are arbitrary complex constants. When it has two different 
stationary points, by shifting the origin and rescaling we can always arrange 
it so that 

W ^) = ^{-^ + ^}- ( 29 ) 

Here we have maintained a phase e tl3 for convenience. The real part and the 
imaginary part of W are 

W R = -- cosP<p R 3 + smfituffa + cos/30^0/ 2 - - sin/30/ 3 

+ cos (3(j) R — sin /3(f)j (30) 

Wi = ~ sin/30 R 3 - cos (3<f)R 2 4>i + sin/30^0/ 2 + ^ cos/50/ 3 

+ sin/30^ + cos/30/. (31) 

In this case stationary points of the superpotential are tp^ = + and 
(f(2) = —-^5 while W(±-^) = ±|e l/3 . Therefore topological charge (19) for 
(2,l)-soliton is 

3(2,1) = +\e iP - (32) 
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(a) (b) 



Figure 1: (a) Gradient curves of W R for (3 — |. (b) Gradient curves of W R 
for 6 = p. 



Comparing this with (22) we have 

e ie = e ip . (33) 



Thus Eqs.(27) become 



and (2,l)-soliton solution is 



d x (f) R = -4> 2 R + 4>] + i 

d^j = 20^0, (34) 



*)r = tanh(x — x ) 

ij = 0. (35) 



Now comparing [Fig. 1(a)] and [Fig. 1(b)], ip^ is not connected with </?( 2 ) 
by any gradient curves of Wr. On the other hand it is connected by a gradient 
curve of Wr. This means that Eq.(25) adjusts gradient curves to connect 
these points by adding Wr and Wj with the ratio of 6. As well the gradient 
curve which connects these points does not depend on 0/. It corresponds to 
4>i = in soliton solution (35). 
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4 Mass Combination Rule 



This section we consider the case where two scalar fields interact with N=l 
super symmetry. The action, including fermions, is 



(36) 



where superpotential H is a real function of fa and 02- The subscript % on 
H(fa, fa) means the derivative with respect to fa 



Hi(fa,fa) 



dH(fa,fa) 
dfa 



(37) 



In this case, minimum points of the scalar potential V = Yji\(Hi) 2 are 
scattered in the (fa, 02)-plane. These points are solutions of equations 



i?i = 
H 2 = 0. 



(38) 



Suppose that there are m solutions 



0(1) = (01(1), 02(1)) 
0(2) = (01(2), 02(2)) 



(39) 



0(m) = (01(m), 02(m))- 

Whenever m > 2, equations of motion have soliton solutions. Boundary 
conditions of (i, j)-soliton are 



\<P1,<P2) = Wl(i),(P2(i)) 
(fa, fa) = (01(j),02(i)) 

The Hamiltonian for the static solution is 



for 
for 



x 



-oo 



x — > +oo. 



n 



J 



dx 



^MfaTH^^dxH 
,i=i 1 



(40) 



(41) 
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(a) (b) (c) 

Figure 2: Gradient curves near the (a) bottom, (b) top and (c) saddle point. 



We derive the Bogomol'nyi mass bound similar to (8), where is 

T (i,j) = #(01(j)> <t>2(j)) ~ #(<£l(i), 2 (i))- (42) 

On this bound, Ti is saturated by solutions of equations 

d x cj)i = ±Hx 

d x <f> 2 = ±H 2 . (43) 

They are similar to (27). Thus again, the integral curve coincides with a 
gradient curve of H. 

However, H is not a harmonic function. Stationary points of H can be 
classified into three types. They may be bottom (minimum), top (maximum) 
or saddle. A bottom point is the source of gradient curves [Fig. 2 (a)] and a 
Top point is the sink [Fig.2(b)]. Namely infinite number of gradient curves 
start from the bottom point, and flow into the top point. Since the soliton 
integral curve coincides with a gradient curve, if the superpotential has a 
bottom point and a top point, equations of motion may have continuous 
infinity of soliton solutions. 

On the other hand, at the saddle point, only two gradient curves can start 
from it and only two gradient curves come into it [Fig.2(c)]. Thus there are 
at most four soliton solutions which pertain to a saddle point on one of the 
boundary conditions. In the previous section, since all stationary points of 
Wr are saddle, there are only finite number of soliton solutions in the N=2 
super symmetry. 
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Now we consider the cubic superpotential case. The general form is 

H = 2 A Hk<t>i<t>j<t>k + B ijMj + Yl C i<i>i ( 44 ) 
ijk ij i 

where A^, Bij, Ci are real coefficients, and symmetric with respect to the 
interchange of indices. Generally in this case, Eqs.(38) have four solutions. 
Among those four points, we can fix two points at (0i,02) = (±1,0). For 
simplicity we further fix remaining two points at (0i,02) = (0, ±p), where 
p > 0. Then H has two parameters: 

1 1 

H = -g0i 3 + a 2 0i 2 02 - ai0i02 2 + ga 2 ai0 2 3 + 4>i - a 2 2 (45) 

where a 2 is an arbitrary constant and a\ = ^ > 0. Now let us call (+1,0), 
(-1,0), (0,+p), (0, -p) as 0(1), (2) , 0(3), 0(4) respectively. 

To investigate the aspect near these points, we consider the eigenvalues 
of the Hessian matrix 



22 



at each stationary point, where 

d 2 H 
50i 
5 2 i^ 

Hn = ki , = 2a 2 0i - 2ai0 2 
#22 = 77—9 = -2ai0i + 2a 2 ai0 2 . 

502 

Of two eigenvalues, if both are positive, it is a bottom point. If both are 
negative, it is a top point. And if one is positive and the other is negative, 
it is a saddle point. See below. 

Ai < 0, A 2 < • • • Top 
Ai > 0, A 2 > • • • Bottom 
AiA 2 < • • • Saddle 

Hereafter we denote the bottom point by B, the top point by T and the 
saddle point by S. 
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Eigenvalues of matrix (46) at 0(i) , (f)(2) , 0(3) , 0(4) are respectively 

A = -( ai + 1) ± v /(a 1 + l) 2 -4(a 1 -a 2 2 ) (47) 



A = (ai + l)± v /(ai + l) 2 -4(a 1 -a 2 2 ) (48) 



A = Q2(v / al+ — !=) ± 



a 2 2 (v^T + -— ) 2 + 4(ai - a 2 2 ) (49) 



/ a v V 1 1 / 



A = -a 2 (Vai + -— ) ± 



a 2 2 (vV + ^) 2 + 4( fll - a 2 2 ). (50) 



The feature of these eigenvalues are dictated by a\ — a 2 2 . Since a± = ^ is 
fixed, a 2 decides the feature of stationary points of H. The classification of 
the types of (0 ( i } , (2) , (3) , (4 )) is (Fig.3) 

(S,S,T,B) for a 2 <-- (51) 

V 

(T,B,S,S) for --<a 2 <- (52) 

p p 

(S,S,B,T) for a 2 > i. (53) 

Values of superpotential i7 at each point are respectively H((f>^) = +| , 
H{4>(2)) = -§ , #(0(3)) = -§P«2 , #(0(4)) = +§P«2- We can check if 
takes the largest value at the top, the smallest value at the bottom and the 
intermediate values at the saddle points. 

Watching (Fig. 4) we find that H have only two gradient curves connecting 
B and S, only two connecting S and T and infinite number of gradient curves 
connecting B and T. On the other hand it does not have the gradient curve 
connecting S and S. 

As an example, let us examine closely case — ^ < a 2 < ^ [Fig.4(b)]. In this 
case, Eqs.(43) have following soliton solutions: (4,l)-soliton, (3,l)-soliton, 
(2,4)-soliton and (2,3)-soliton, and infinite number of (2,l)-solitons. But 
there is no (3,4)-soliton. Further (2,l)-soliton solutions make one parameter 
family with a common mass. It may be considered bound states either of 
(2,3)-soliton and (3,l)-soliton, or (2,4)-soliton and (4,l)-soliton. And masses 
of these solitons satisfy an important relation 

M 21 = M 23 + M 31 = M 24 + M 4l . (54) 
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ai 



ai=(a 2 ) 



(T,B,S,S) 



(S,S,T,B) \ 


/ (S,S,B,T) 




32 



(S,S,-,-) 



Figure 3: Classification of the type of stability points (4>(i), <fi( 2 ), 0(3), 0(4))- 



This reminds us of the combination rule of Ritz and implies a relation of 
marginal stability for (l,2)-solitons. In our special case, masses of these 
solutions are M 41 = M 23 = |(1 — pa 2 ) , M 31 = M 24 = |(1 + pa 2 ) and 
M 21 = §. In the case of a 2 < -- or a 2 > -, the relation of soliton solutions 

is similar to the case of — - < a 2 < -. 

p p 

Case a± = 1 and a 2 2 ^ 1 is special. In this case, stationary points of H 
are at (0i, fa) = (±1, 0), (0, ±1). We can solve Eqs.(43) as follows. In terms 
of new variables 

1 



+ fa} 



7=2^ 



(55) 

^2}, (56) 
superpotential (45) becomes 

H = V2 {(1 - a 2 ){- l -4> + * + + (1 + a 2 ){-^J + , (57) 

which is the sum of <p + part and </>_ part. Equation (43) reduces to two 
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separate equations whose solution is 

+ = ±-^tanh((l-a 2 )(x-:4)) (58) 
0_ = ±-^=tanh((l +a 2 )(x-Xo)) , (59) 

where Xq or Xq is the center of each soliton. The parameter of the soliton 
solution curve which connect B and T is ^(xq—Xq). And \{xq +Xq ) indicates 
the center of 0+ soliton and 0_ soliton. 

We have obtained two non- interacting solitons posed at rest. Each soliton 
belongs to different degree of freedom: + and 0_. Its obvious generalization 
to iV scalars = (0i, 02, • • • , 0jv) is 

# (0 = E (3^ - • ( 6 °) 

Models which we considered in this section have many properties in common 
with this trivial case. We do not know, however, general criterion when this 
sort of disentanglement of superpotential occurs. 



5 Promotion from N=l to N=2 Theory 

In this section, we will investigate when and how the promotion from N=l to 
N=2 takes place in two coupled scalar fields. As the superpotential in N=2 
theory is holomorphic, we are hinted to force our N=l superpotential to be 
harmonic: 

( d 2 d 2 \ 

Then we get a\ = —1, and 02 is arbitrary. In our previous context, this is 
put that p must be ±i It is curious that we obtained harmonic H, not by 
pushing 0(3) and 0(4) to infinity, but by making them imaginary. 

Now, in order to understand what happens at a\ = —1, let us assume 
a\ < 0, but not necessary a\ = —1. If 0i and 02 regard as 4>r and 0/ 
respectively, we see that a 2 correspond to tan/3 in (30) or —^-3 in (31). 
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Therefore let us replace a 2 by tan (3, and rename H cos (3 as H. 

H = — -cos/?0i + sin/3^1 02 — ai cos/30102 +-aisin/30 2 
o o 

+ COS/901 - sin/302 (62) 

Then all stationary points are saddle points (Fig. 3) from (47) and (48). 

When /3 ^ and (3 ^ n, 0( X ) is not connected with 0( 2 ) by any gradient 
curves [Fig. 5(a)]. All gradient curves that start from or come into stationary 
points extend to infinity. Therefore Eqs.(43) do not have soliton solution. 

When (3 = 0, H becomes [Fig.5(b)] 

^ = -^0i 3 -a 1 1 2 2 + 0i. (63) 

One of the gradient curve connects 0(i) and 0( 2 ). Thus there is a soliton 
solution. When (3 — 7r, the result is the same as this. 

However when Oi = —1, the situation is different. In this case H(<j>i, 2 ) 
becomes harmonic function, and it has a conjugate function ft,(0 l5 2 ) de- 
fined by Cauchy-Riemann equations. Therefore we can define a new phase a 
between these functions. Subsequently N=l supersymmetry is promoted to 
N=2. If i/(0i,02) has the same form as Wr (30) when (0_r, 0/) is replaced 
by (0i, 02), its conjugate function /i(0i, 02) is given by Wj (31). 

Then we can rewrite the scalar potential in (36) 

£ \(H) 2 = £ \{Hi cos a + hi sin a) 2 (64) 
i=i 1 i=i 1 

and (41) become 

r 2 1 

Ti — / dx[ —{d x 4>j — (Hi cos a + hj sin a)} 2 
J i=i 2 

-l-c^i? cosa + /isina) ]. (65) 

The mass bound is 

Ti > T(jj) cos a + t(j j) sin a, (66) 

where is 

= h (<t>(j)) ~ M </>«)• (67) 
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This bound is saturated by solutions of equations 

d 

d x (f)i = -^—-(H cosa + hsma) (68) 
d 

d x 4>2 = tt— (Hcosa + ft, sin a) (69) 

0<p2 

When a = (3, H cos a + h sin a is the same as H with /3 — 0. When ct = /3+7T, 
it is the same as if with /3 — ir. That is, by adjusting a, we can always put 
/3 = or /3 = 7r. Thus when ai = —1, equations of motion have always 
soliton solutions. These discussions are consistent with Sec. 3. This means 
that N=l supersymmetry is promoted to N=2 when a\ = —1. 

We summarize ranges where equations of motion have soliton solution. 
They are a± > 0, a± = —1, and a\ < with 02 = 0. 

Let us notice that we can define the conjugate function of H for arbitrary 
ai, although they do not satisfy Cauchy-Riemann equations. It is derived 
from H (62) by replacing f3 to f3 — |. 

h — — sin/30! — cos/3^1 02 — «i sin/30i02 aiCos/30 2 

3 3 

+ sin /30i + cos /30 2 (70) 

We can introduce the phase a between (62) and this, and Hcosa + hsina 
is same as H with (3 = (63). This construction is similar to the harmonic 
case. But since they do not satisfy Cauchy-Riemann equations (26), we can 
not deform the scalar potential like (64), and we can not put (3 = or /3 = n 
by a. 



6 Summary 

We considered the soliton in two dimensional N=l supersymmetric theory in 
which two fields couple through quasi-superpotential. In this case, masses of 
soliton solutions satisfy combination rule 

M12 + M 23 = Mi 3 . 

We demonstrated this result along with special two parameter family 
of superpotential (45). However, this combination rule of soliton masses is 
considered a general feature of non-harmonic superpotential. 
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This implies marginal stability of solitons, and suggests the absence of 
intersoliton force between some solitons in this model. To check it by several 
ways will be our future work. 
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